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Abstract 

This report deals with instrumentation and computer 
progr dumnj concepts that have been developed for ultrasonic 
materials characterization. Methods that facilitate velocity 
and attenuation measurements are described. The apparatus 
described is based on a broadband, oufferea contact prone 
using a pulse-echo approach for simultaneously measuring 
velocity and attenuation. Instrumentation, specimen 
condition, and signal acquisition and acceptance criteria are 
discussed. Typical results with some representative materials 
are presented. 


1. INTRODUCTION 
1. 1 BACKGROUND 

There is high interest in the potential 
uses of ultrasonics for nondestructive 
characterization of materials. Physical 
bases and a rationale for quantitative 
ultrasonic evaluation of mechanical 
properties of materials ere recently 
reviewed <*>. It was noted that even 
where conventional nondestructive 
evaluation, NDE, techniques have shown 
that a part is free of overt defects, 
advanced NDE techniques should be 
available to confirm critical material 
properties assumed in the part's design. 

Aaoag the available NDE techniques 
(including radiometry, eddy currents, 
■icrowaves, etc.) ultrasonics is now 
proainent as a highly versatile means 
for mechanical property character- 
ization. Mechanical properties that can 
be ultrasonically evaluated include 
elastic moduli, tensile strength, yield 
strength, ductility, hardness, and 
fracture toughness. 

Ultrasonic measurements of mechanical 
properties depend on making accurate and 
precise velocity and/or attenuation 
measurements. Examples where both kinds 
of measurements are used include 
determination of hardness and tensile 


strength *•» and correlations of yield 
strength and fracture toughness <*>. 

1.2 ORGANIZATION 

This ceport describes instrumentation 
and programming concepts that facilitate 
velocity and attenuation measurements in 
ultrasouic NDE of material properties. 
Practical factors essential to making 
ultrasonic measurements for quantitative 
material property characterization are 
discussed. Operational details are 
given as a guide for further system and 
instrumentation developments. 

The waveform aquisition and processing 
concepts described generally relate to 
laboratory material samples. However, 
the instrumentation and programming 
guidelines are applicable to a variety 
of actual structural parts if specified 
size, geometry, surface guality, and 
morphological conditions are satisfied. 

This report has four major sections: 
INSTRUBENTATION, SIGNAL ACQUISITION, 
SIGNAL PROCESSING, and REPRESENTATIVE 
RESULTS. 
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2. I NSTRUHENT ATION 

2.1 OVERALL STSXEB 

An ultrasonic waveform processing systen 
based on a high Ereguency digitizer and 
high speed minicomputer was asseabled 
froa coaaercially available components. 
At this writing, the systeo appears to 
be unique in its capaoilities for per- 
forming high frequency ultrasonic signal 
ana lyses. 

A contact probe, aultiple pulse echo 
approach is used for the velocity and 
attenuation aeasureaents. Simplicity of 
application and versatility are afforded 
in combining both the pulse emitter and 
signal receiver in a single ultrasonic 
probe that requires access to only one 
surface of a sample material. The 
signal acquisition and processing 
arrangement is shown in figs. 1 to 3. 
Tfca principal system components are: 

(1) Buffered ultrasonic transducer 

(2) Broadband pulser-receiver 

(3) Programmable waveform digitizer 

(4) Interactive graphics terminal 

(5) Minicomputer and software system 

The above major components and pertinent 
operational factors are described below. 

2.2 ULTRASONIC TRANSDUCER 

It is advantageous to use a buffered, 
broadband ultrasonic transducer such as 
that illustrated in fig. 1. The buffer 
isolates the transducer crystal from the 
specimen and eliminates adverse loading 
effects that can disturb and distort 
incoming signals. 

A single 50 MHz transducer is the basis 
for data reported herein. It contains a 
car amir piezoelectric element polarized 
to emit longitudinal waves. The buffer 
is a fused quartz circular solid 
cylinder with polished ends. The buffer 
length is 12.7 nm (0.5 in) giving a 
buffer delay fixed at 4.3 nicroseconds. 
The buffer diameter is 16 mm (.63 in) 
and the piezoelectric crystal has a 6. 3 
nn (.25 in) diameter. 

2.3 PULSER-R ECBI VBR 

The ultrasonic pulser- receiver is a 
special manufacturer- modified unit. 
The aodifications consist mainly of a 
sharpened excitation pulse and extended 
high frequency drive in order to achieve 
greater aaplification and wider bandpass 
perforaaace. The receiver amplifier has 
a 75 MBz bandpass. 


The excitation pulse has a 200 volt 
peat, 3 nanosecond risetiae, and 10 
nanosecond duration. These are minimum 
values that can be increased by an 
energy control dial. The instrument 
also includes a damping control for 
inpedance matching the transducer with 
the pulser- receiver and an attenuation 
control for reducing the effects of 
distortion in the amplifier circuit. 

2.4 WAVEFORM DIGITIZER 

Waveform acquisition for ultrasonic 
velocity and attenuation measurements is 
accomplished by means of a high 
frequency waveform digitizer in which 
analog time domain signals are converted 
into digital arrays. The waveform 
digitizer is a state-of-art controller- 
based instrument with a silicon diode 
scan converter matrix. 

High gain and bandwidth are achieved by 
signal capture on the diode target 
matrix. Signal capture is done at high 
speed while subsequent digitization 
proceeds more slowly at computer 
compatible data transfer rates. Signals 
containing frequencies from 10 Hz to 1 
GHz are readily digitized by the scan 
converter. 

Signal and synchronization inputs to the 
digitizer are made through modular 
oscilloscope plug-ins consisting of a 
dual channel vertical amplifier and a 
time base with delayed and delaying 
sweep provisions. The vertical module 
operates at a 300 MHz bandwidth. The 
time base has a sweep rate of 500 
picoseconds per division. This allows 
signal capture within a 5 nanosecond 
time frame. An x-y video oscilloscope 
monitor is used to view signals in 
real-time as they are written on the 
diode target matrix. 

The digitizer produces 512 element 
waveform arrays. After scan conversion, 
raw waveform arrays are retrieved under 
program control from a self contained 
(4096-word) memory. Each of the 512 
waveform array elements is associated 
with a horizontal time value and 
vertical voltage deflection representing 
the time domain signal. 

2.5 GRAPHICS TRRHIN AT. 

In addition to the usual programming and 
data display operations, the graphics 
terminal is a vital link for interactive 
signal processing. Raw or processed 
waveforms, program generated graphics, 
and other data displays are called up by 
direct operator intervention or under 
program control. This enables operator 
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interaction with intermediate as well as 
final results, when necessary. 

The systen includes a hard copy unit for 
reproducing program listings, data, and 
graphics from the terminal screen. 
Documents are deliverable at about 10 
second intervals. Example of figures 
produced by the graphics package appear 
throughout this report. 

2.6 CONTROLLER AND SOFTWARE 

The controller consists of a high-speed 
minicomputer with a 64,000 word memory 
serviced by two 1.2 million word disks. 
The controller has a real-time operating 
system, memory management, universal 
bootstrap, automatic self diagnotics, 
and floating point processor. Of the 
64,000 word memory, about 20,000 words 
are occupied by the operating system, 
about 12,000 words are reserved for 
programs and processing, and the 
remaining 32,000 words are available as 
"virtual" memory for expanded high speed 
data and array processing. 

The computer resident software system 
incorporates a high level language 
(modified BASIC) < ] > specifically 
evolved for the waveform digitizer by 
the manufacturer's software development 
group. The software includes modules 
for instrument control, peripheral 
support, graphics, waveform and array 
processing. The waveform and array 
processing nodules include suoroutines 
for convolution, cross correlation, auto 
correlation, signal mean, maximum and 
minimum, zero cross, root mean square, 
differentiation, integration, and "fast" 
Fourier transforms "FFT" <*>. 

3. SIGNAL ACQUISTION 

3.1 GENERAL CONDITIONS 

This section describes factors that 
govern the acquisition of viable 
signals. Guidelines are suggested for 
material preparation and probe coupling. 
In addition, procedures for signal 
adjustment, definition, and acceptance 
relative to making accurate attenuation 
and velocity measurements are provided. 

The signal definition and acceptance 
criteria are predicated on materials 
that are linear elastic, nondispersi ve, 
and essentially isotropic. Haterials 
are assumed to have a uniform, poly- 
crystalline, fine grained microstructure 
with the mean grain size being much less 
than the mean wavelength of the probe 
signal. This assures operation in the 
scatter attenuation regime which is 
sssumed for the attenuation coefficient 
measurements described herein. 


3.2 MATERIAL SAMPLE 

To accomplish the quantitative 
attenuation and velocity measurements 
needed to measure material properties 
ultrasonically, the material sample must 
satisfy some specific size, shape, and 
surface quality constraints. 
Acq uisition of optimal signals requires: 

(1) Clean, mach ine- fi nished surfaces 

(2) Flat, smooth, parallel surfaces 

(3) Precisely machined thicknesses 

(4) Large area relative to probe size 

(5) No internal defects or anomalies 

Flat, ground or machined surfaces assure 
good probe coupling. For absolute 
velocity and attenuation measurements 
ultrasound pulses should be introduced 
at normal incidence. Signal trans- 
mission should be between parallel 
surfaces. These constraints reduce the 
deleterious effects of is am divergence 
and other geometry related factors. 

The distance between the parallel 
surfaces used in multiple pulse echo 
measurements should be precise and 
relatively small. If the material 
thickness is too great, attenuation 
losses will tend to diminish the chance 
of recovering strong echoes, 13 1 and B2. 
Conversely, material sample thickness 
should be great enough to give a good 
basis for velocity computations and to 
include a fair sampling of micro- 
structure. 

The region to be probed should be large 
relative to the probe beam diameter and 
free of facets or boundaries that will 
cause extraneous reflections. The 
region being probed should also have an 
ample area to permit probe shifting to 
alternative positions to get good 
lateral sampling of the material. 

Ooviously, other NDE methods should be 
used to verify the assumed absence of 
overt defects, internal anomalies, or 
material gradients that might make the 
sample unrepresentative of virgin 
material. 

[Specimens with plane parallel surfaces 
with smooth finishes have been used. 
Surface roughness was less than 16 rms. 
The specimen surfaces were machine 
ground and subsequently polished with 
600 and 800 grit sandpaper. Typically, 
the lateral area was from 2 to 4 square 
cm. ] 

3.3 BUFFER MEDIUM 

The buffer diameter should as « rule be 
roughly two to three times that of the 
crystal to avoid sidewall reflections. 



The length should be such that at least 
two back surface echoes from the 
specimen, B1 and 62, Ray be readily 
included between successive buffer 
echoes, FI and F2, as in fig. 2. If 
this condition is net, it is easier to 
identity and isolate echoes 01 and B2. 
Otherwise, these echoes can becone 
confuses with or run into other nultiple 
echoes in the time domain. 

Fused silica and quartz are equally 
appropriate buffer materials. These 
materials have a vitreous, amorphous 
microstructure with very low attenuation 
in the frequency range of interest. 
Also, the acoustic impedance ot either 
fusel silica or quartz is appropriate 
for the materials studied. The buffer 
acoustic impedance should be less than 
that of the material to which it is 
coupled. These factors assure ample 
signal amplitude for the echoes returned 
by the material, B1 and B2, in fig. 2. 

[The buffer in the 50 HHz transducer was 
fused quartz. The buffer is usually an 
integral part of commercial probes and 
its length is fixed by transducer design 
requirements. Therefore, the specimen 
thickness was coordinated with buffer 
length to assure two material echoes in 
the buffer delay interval. For example, 
for the materials studied the specimen 
thickness ranged from 2 to 3 mm. ] 

3.4 COUPLING MEDIUM 

Coupling media such as glycerine, light 
oils or commercial gel couplan+s will 
perform well if the previously noted 
material surface conditions are net. 
The couplant should readily wet the 
material surface. 

The couplant thickness should be a 
minimum and highly uniform. Otherwise, 
multiple reflections and interference 
effects within the couplant tend to 
degrade any high frequency signal 
components being transmitted. This is a 
serious undesirable effect in the case 
of attenuation measurements that rely on 
spectrum analysis of high frequency 
signals. The couplant thickness should 
be reduced to a minimum by squishing. 
The application of appreciable pressure 
is needed for this. 

[For the glycerine couplant used in this 
study a squishing force of from 44 to 88 
N (10 to 20 lbs) was found appropriate. 
This was equivalent to from 220 to 440 
kPa (30 to 60 psi). Specimens were 
supported atop an instrumented load cell 
tor these measurements. ] 


3.5 PULSER-RECEI VER ADJUSTMENTS 

As with all instruments of this type, 
impedance matchiug of the pulsep- 
receiver circuitry and transducer is a 
vital factor. The interconnecting 
cable's length and impedance should be 
selected to produce optimum signals. To 
further optimize performance of the 
transducer- instrument combination, fine 
control of impedance matching is 
afforded by a dampinj control in the 
cable circuit. Under proper conditions 
this control will reduce or virtually 
eliminate "ringing" or reverberations by 
altering impedance loading between the 
transducer and pulser- receiver. 

Although the resultant signal may be 
weaker, it is preferable to use minimum 
pulse energies and hence lower pulse 
widths with high frequency transducers 
to get the greatest possible bandwidths. 
A pulse energy triming control in the 
pulser circuit should be available for 
this adjustment. 

An attenuation adjustment in the 
receiver amplifier circuit should be 
increased as required to reduce signal 
distortions. These may arise due to 
amplification nonlinearities. It may be 
necessary to reduce signal strength of 
echoes U1 and B2 while attempting to 
optimize the waveform. 

3.6 SIliNAL RESOLUTION 

To optimize the digital version of the 
time domain waveform, back echoes 61 and 
B2 should be separately digitized into 
their corresponding arrays. Figure 4(a) 
shows a composite waveform sythesized 
into a 1500 element array from the two 
separate 512 element arrays in figs. 
4(b) and 4(c). The purpose of the 
composite waveform is solely to 
illustrate the relatively long tine 
delay between successive echoes, Bl and 
B2. That is, the delay between echoes 
is very long relative to the duration of 
oscilla tions. This is typical for high 
frequency signals in the material 
samples studied. 

Any digitiza ion of ultrasonic signals 
that compresses two echoes into a 512 
element array will greatly reduce fine 
structure definition of the waveform of 
either echo. It is desirable to reduce 
the data sampling interval to the lowest 
possible value for each echo to be 
sampled. Accordingly, the echo waveform 
should "fill the screen" (i.e. , scan 
converter matrix) in order to get the 
best signal resolution. 
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[Typically, the tine domain digital 
saipling interval was from 2 to 0.4 ns 
(between adjacent tiae points of the 
waveform array). This required sweep 
rates of froa 1000 ns to 200 ns per 
wavefora. ] 

3.7 SIGNAL ACCEPTANCE 

A programmed pattern recognition routine 
or simple visual inspection should be 
used to reject unacceptaole waveforms. 
Preliminary critical exaaination of the 
signals should be done each tine the 
probe is recoupled to the material 
sample. 

Certainly, the signals should be free of 
random spikes or any significant noise. 
The tiae and frequency domains of an 
ideal wavefora is illustrated in fig. 
5(a). AcceptaDle waveforms are shown in 
figs. 5(b) and 5(c). Unacceptable 
waveforms appear in figs. 5(d) through 
5(f). 

Within the tine domain, an acceptable 
wavefora is fairly symmetrical and has a 
M W M shape. The two successive negative 
oscillations should have nearly equal 
minima. In the polar version of the 
frequenc domain, the waveform spectrum 
should exhibit a nearly Gaussian shape. 
The spectrum should exhibit no prominent 
multiple peaks, valleys, or distortions. 

Acceptable waveforms will occur if the 
previously mentioned conditions and 
specimen constraints are satisfied. 
Unacceptable or much less than ideal 
waveforms nay occur because of: 

(1) Defective probe/crystal/buffer 

(2) Poor impedance matching 

(3) Binging, underdanping 

(4) Poor/nonuniform coupling 

(5) Amplification nonlinearities 

(6) Scan converter nonlinearities 

(7) Material anomalies/defects 

Generally, waveform problems are easily 
remedied by instrumentation adjustments 
and better sample preparation and higher 
coupling forces. Good signals will tend 
to "pop in" and stanilize when optimum 
coupling is achieved. 

4. SIGNAL PBOCESSING 

4.1 OBJECTIVE 

In many instances either velocity or 
attenuation measurements alone suffice 
for evaluating material properties. 
However, both measurements are necessary 
to evaluate fracture toughness and yield 
strength <*>. In this case velocity 
measurements are needed to determine 
acoustic impedances and to augment 


attenuation measurements from which 
attenuation coefficients aru found. The 
variation of the attenuation coefficient 
as a function of ultrasonic frequency is 
the basis for deducing material fracture 
toughness and yield strength. 

The object, therefore, is to illustrate 
programmed waveforir. processing for 
velocity and attenuation measurements 
leading to fracture toughness and yield 
strength correlations. Digital signal 
processing procedures given herein are 
based in part on the analog procedures 
described in ref. 2. 

4.2 INITIALIZATION AND INTERACTION 

Although the waveform acquisition and 
analysis are performed under computer 
control, certain operator actions are 
required to start and oversee the data 
processing. The establishment of probe 
coupling, for example, is done manually. 
Data concerning the particular material 
samples being probed are input manually 
at the terminal keynoard. As the signal 
processing proceeds, the operator is 
prompted for various additional inputs. 
Operator inputs include: 

(1) Specimen identification 

(2) Thickness and density 

(3) Average grain size 

(4) Coupling medium 

(5) Transducer specifications 

(6) Instrument settings 

(7) Probe pcsition information 

(8) Signal acquisition options 

(9) Signal processing options 

Specimen data includes a serial number 
and material or alloy identity. It also 
includes any available previous property 
measurements such as hardness, yield 
strength, toughness, etc. 

Previously determined material sample 
density and thickness values are input 
at the initialization stage. An ASTH 
Stardard grain size determined for each 
material sample is also input. However, 
this can be circumvented for materials 
for which a relation between velocity 
and mean grain size has been previously 
established <*>. In this case, an 
expression for the functional relation 
between velocity and mean grain size may 
be loaded into the program and the grain 
sized is computed from the velocity 
measu cements. 

Transducer specifications include: the 

nominal (center) frequency, wave node 
(transverse or longitudinal) , buffer 
material. Given the buffer material, 
the program assigns an appropriate 
acoustic impedance value for use in 
computing reflection coefficients. The 
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acoustic iapedaace of the couplant is 
also determined from previously stored 
data. 

Key pulser- receiver settings and other 
nonprograaaabie instruaent settings are 
input for documentation purposes. Prone 
positioning information is requested by 
the program when: (a) the probe is 
repositioned, or (b) the probe is kept 
stationary for signal averaging. This 
enables repetitive sampling at a given 
position or multiposition sampling. The 
program will return a tabulation of 
individual position and averaged data 
values. These data include velocity and 
attenuation parameters that can be 
correlated with material properties, as 
explained under BEPBESENTATIVE BESULTS. 

Opon startup, the operator makes 
preliminary judgements of acceptability 
from the real-time signals displayed on 
the x-y monitor. As the waveform 
processing progresses, the operator may 
exercise further options based on 
program generated graphics and data. 
The options include a selection of 
alternatives from a repertoire of 
auxiliary signal processing routines. 

4.3 VELOCITY NEASOBEBENT 

Velocity is calculated from a simulated 
"pulse overlap" routine in which 
corresponding parts of the successive 
echoes, B1 and B2, are compared and 
matched. Using the previously input 
material thickness, velocity is computed 
from the mean tine delay between the 
echoes. 

The method of pulse echo overlap for an 
analog system is described in refs. 2,6. 
In an analog system two succesive 
material echoes are manually overlapped 
on a dual display oscilloscope screen. 
The tine delay reguired to bring key 
parts of the two echoes into coincidence 
is the basis for computing velocity. 

The overlapping and velocity computation 
are performed automatically by the 
waveform processing system using digital 
waveform and array processing routines. 
Array processing routines find specific 
tine elements associated with criterion 
points of the waveform. Criterion 
points inclnde the peak values of the 
first and second major oscillation and 
the xero crossing point between then. 
The criterion points are located by an 
array search process. The mean time 
delay between successive echoes, B1 and 
B2, is found by determining the time 
delay between corresponding criterion 
points. This is accomplished by 
reference to a continuous pulse train, 
ns described nest. 


4.4 REFEBEHCE PULSE TBAIN 

A precision, pulse generator is used to 
produce a continuous pulse train that 
alternates with the ultrasonic signals 
on the x-y monitor and scan converter 
matrix, see fig. 6. The ultrasonic 
signal is fed to Channel 1 while the 
pulse train is fed to Channel 2 of the 
dual trace vertical amplifier. 

Since high signal resolution demands 
that only one ultrasonic echo at a time 
be digitized, the reference pulse train 
is a necessary link between the echoes. 
The pulse train serves as a time marker 
by means of which any arbitrary delay 
adjustment needed to bring either B1 or 
B2 into the (scan converter) screen area 
may be determined. The pulse period and 
phase are adjusted so that a specific 
pulse "X" will appear with B1 in its 
tine frame, while a specific successive 
pulse ”Y" will appear with B2 in its 
tine frame. 

The pulse train period is read into the 
waveform processing program. Once the 
two succesive echoes have each in turn 
been positioned on the (scan converter) 
screen by the time base delay control, 
the delay is computed from, 

(TD) = (TP) - (TX) ♦ (TY) 

where, (TD) is the time delay reguired 
to overlap the echoes, (TP) is the pulse 
train period, (TX) is the signed tine 
interval between pulse "X" and the mean 
criterion point of Bl, and (TY) is the 
signed tine interval between pulse "Y" 
and the mean criterion point of B2. The 
algebraic signs of (TX) and (TY) are 
positive when their corresponding 
pulses, X and Y, respectively, lead the 
echoes and negative if they lag. This 
is is illustrated in fig. 6. The 
program then computes velocity (VL) 
based on the delay from, 

(VL) * 2 (ST)/ (TD) 

where, (ST) is sample thickness. 

[Throughout this report the upper case 
letters enclosed within parentheses 
represent symbols that always denote 
specific quantities: frequency (P) , 
reflection coefficient (BC) , attenuation 
coefficient (ALPHA) , slope paraaatwr 
nunber 1 (BETA 1 ) , and so forth.] 

4.5 ATTEB UATION HEASOBEHENT 

Bhile velocity reaains fairly constant 
relative to frequency (usually varying 
lass than one percent) over a wide ramge 
of frequencies, attenuation will 
generally vary by several orders of 


nagnitude for the sue frequency range. 
Frequency spectrua analysis is the basis 
for deteraining the variation of 
attenuation as a function of frequency. 
The frequency spectrua analysis aethod 
for an analog systea is described in 
refs. 2,6. 

The objective is to find attenuation 
coefficient (ALPHA) as a function of 
frequency using the equation, 

(ALPHA) * (LOGN[ (BC)/(BB) ]) / [2(ST)] 

where, (BC) is reflection coefficient, 

(BB) is the aaplitude ratio B2/B1 at 
specific frequencies over a range of 
frequencies, and LOG N denotes natural 
logaritha. Both (BB) and (ALPHA) are, 
of course, functions of frequency while 

(BC) is essentially independent of 
frequency. The above equation applies 
only if the opposite side of the 
aaterial saaple is free, that is, 
exposed to a gaseous aediua like air. 

To find (ALPHA) as a function of 
frequency (F) , it is first necessary to 
deteraine (BC) and a functional relation 
between (BB) and (F) . This is done by 
aeans of frequency spectrua analysis and 
the aaplitude ratio curve, ABC, as 
described next. 

4.6 FREQUENCY SPBCTBUH ANALYSIS 

The ultrasonic frequency spectra are 
obtained by the discrete fast Fourier 
transforaation, FFT, of tiae doaain 
signals. Echoes B1 and B2 are 
individually analyzed into frequency 
doaain arrays. These arrays yield the 
raw aaplitude ratios BVB2 over a range 
of frequencies of interest. The ratios 
are conputed at 1 MHz intervals froa 0 
to 100 HHz. 

The FFT is perforaed on two versions of 
each of the echo waveforas. First, the 
FFT's of the natural waveforas of B1 and 
B2 are created. Second, the FFT is 
perforaed on truncated versions of B1 
and B2. The truncation consists of 
setting ainor trailing oscillations in 
the wavefora arrays equal to zero. 
Thus, in the truncated version only the 
first five aajor oscillations of the 
classic "N" pattern are retained. The 
natural and truncated versions of the 
signal are congruent except for the 
excluded tail oscillations. 

In the frequency doaain, the aaplitude 
spectra of the noraal and trucated 
versions are approxinately congruent, if 
the original signals aeet the acceptance 
criteria given previously. This is 
illustrated in fig. 7(a) where the solid 
corves are for the natural waveforns and 


the dotted curves are for the truncated 
versions of B1 and B2, respectively. 

Not all the frequency doaain data is 
usable. A zoning procedure is required 
to define a frequency range iron which 
valid attenuation data aay be extracted. 
Zone linits are set by aeans of the 
first and second derivatives of the 
frequency spectrua of the weakest echo, 
B2. The differentiation is perforaed on 
the truncated version since it gives 
aore definitive results. 

The aaxiaua of the first derivative, 
fig. 7(b), defines the lower frequency 
liait of the zone. This aaxiaua occurs 
oelow the center frequency of B2. The 
second peak of the second derivative 
defines the zone high liait. This peak 
occurs above the center frequency of b2. 
Zone linits defined in this way are 
indicated by vertical narks. It is seen 
that the linits correspond to critical 
infelctions in the spectra. 

4.7 AMPLITUDE BATIO CURVE, ABC 

The data points in fig. 8 are a plot of 
spectral aaplitude raw ratios 82/Bi 
versus frequency. The solid curve 
fitted to the raw data is the aaplitude 
ratio curve, ABC. The ABC represents 
the functional relation between (BB) and 
(F). It is evident that the ARC is 
fitted to only the zoned portion of the 
raw data. The vertical Barkers indicate 
the valid zone defined according to the 
previous criteria. 

Excluded raw data to the left of the 
valid zone are within the diffraction 
regiae <»». This regiae typically 
extends to aoout 15 or 20 HHz. It is 
possible to adjust the raw data by aeans 
of diffraction corrections. However, 
the need for any diffraction corrections 
is avoided by having abundant higher 
frequency data <•>. 

External raw data to the right of the 
valid zone are discarded chiefly because 
they fall outside the effective 
bandwidth sensitivity range of the 
instr uaentation (that is, the pulser- 
receiver- probe coabination) . 

Using only the raw data froa the valid 
zone, a functional relation between (BB) 
and (F) is found by an iterative 
anltiple regression routine. This 
involves three key paraaeters: the 
coefficient (C) , the exponent (H) , and 
reflection coefficient (BC) in, 

(LOG»[ (BC) / (BB) ]| / [ 2(ST) ] « (C) (F) 
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The previous equation is based on the 
relation between attenuation coefficient 
(ALPHA) and frequency (P) given by, 

(ALPHA) * (C)(P)‘*» 

This expression for (ALPHA) is valid 
under the conditions that scatter 
attenuation predominates and that (C) 
and (H) can be taken as constants over 

the frequency range of interest <*>. 

These conditions Apply under the 
■aterial constraints stated in GENEBAL 
CONDITIONS. 

4.8 B EPLECTI ON C0EPP1CIENT 

The key to fitting an ABC to the zoned 
data depends on finding the reflection 
coefficient (HC) for the buffer- sanple 
interface. Because of vitiations in the 
degree of coupling, the reflection 
coefficient (BC) is not a fixed or known 
quantity. If the buffer and naterial 
sanple interfaced perfectly, as under 
very high pressure with no intervening 
couplant, the reflection coefficient 
would be, 

(BCB) * [ (ZS)-(ZB) ] / [ (ZS ) ♦ (ZB) ] 

where, (ZS) is the sanple naterial 
acoustic inpedance and (ZB) is the 
buffer acoustic inpedance. 

Hhen a coupling nediun is introduced to 
get acoutic continuity, (BC) will differ 
fron (BCB). Generally, (BC)a(BCB), 
given that (ZS)>(ZB), as required by the 
conditions nentioned earlier under 
BDPPEB MEDIUM. (BCB) is thus a lower 
bound on (BC) . Also, (HC) < (BCC) , where, 

(BCC) « [ (ZS)-(ZC) ] / [ (ZS) ♦ (ZC) ] 

(ZC) is the couplant acoustic inpedance. 
The quantity (BCC) is thus an upper 
bound on (BC) . 

He define (BC)MBCV) where (BCV) is the 
"valid" reflection coefficient. Then, 
(RCB) S (BCV ) < (BCC) . (BCV) is deternined 
by iterative regression analysis of the 
zoned raw data for B2/B1 vs. frequency. 

The iterative regression routine assigns 
a starting value to (BCV) equal to 
(BCC). A two paraneter multiple 
regression routine finds a succession of 
interia values for (C) and (N) . A 
standard deviation is conputed by 
conparing array values of the zoned raw 
data with interia ABC's conputed fron 
the current values of (H) , (C), and 

(BCV). The regression proceeds for a 
series of decreaented (BCV) values. Por 
each iteration, a test is performed on 
the current standard deviation. This 
iterative process ends shen: (a) the 


current standard deviation value is a 
nininun, or (b) (BCV) S (BCB). The 
process repeats until the value for 
(BCV) is deternined to a precision of at 
least four significant figures. An 
exanple of the iterative curve fitting 
is shown in fig. 9. 

Two tests are nade to deteraine the 
acceptability of the data at this stage: 
First, if (BCV) < (BCB), all waveforns and 
associated data are rejected. Second, 
rejection also occurs if the coefficient 
of deteraination or goodness of fit 
(PIT) is less than a specified value 
(e.g., <0.99 6 ). In either case new 

waveforns aust be acquired and the 
processing reinitiated. 

5. BEPBESENTATIVE BESOLTS 

5.1 ATTENUATION VS. FREQUEMCI CURVE, AFC 

Upon completion of the the ABC fitting 
process, the prograa returns values for 
(BCV), (M) , and (C) . These quantities 
are displayed and docunented as shown ny 
the exanples in fig. 10. The solid line 
superiaposed on the raw data in fig. 10 
is the ABC or (BB) versus (P) found 
fron, 

LOGN(BB) = LOGN(BCV) - 2 (ST) (C) (P) 

Figure 10 shows results obtained with 
both versions of the wavefora: the 

natural and the truncated. Differences 
that are evident in the corresponding 
raw ratios B2/B1 are explained under 
DISCUSSION. 

Establishnent of the ABC yields the 
paraneters needed to define the 
functional relation sought between the 
attenuation coefficient and frequency. 
Sone typical plots for the attenuation 
vs. frequency curve, APC, appear in fig. 
11. The solid curve superiaposed on the 
raw data is based on the values for (H) , 
(C) , and the equation, 

(ALPHA) = (C)(F)<*> 

The APC characterizes the naterial 
sanple aicrostructure and, therefore, 
certain naterial properties. Further 
analysis is needed to use infornatioa 
contained in the APC. 

5.2 HATEBIAL COBBELATIONS 

Correlations will exist with properties 
that are governed by the sane 
aicrostructural factors that govern the 
velocity and attenuation paraneters. 
Por exanple, it has been found that 
certain naterial properties correlate 
closely with velocity and attenuation 
factors that can be derived fron the APC 
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parameters. The significance of the 
previous signal processing results is 
explained in terns of correlations with 
fracture toughness and field strength. 
The approach is to check the results 
against independently seasured aaterial 
properties obtained previously for a set 
of standard specinens. The standard 
rpeciaens in the present case consist of 
250-grade aaraging steel. Prior 

aechanical tests established fracture 
toughness and yield strength for these 
specimens. 

5.3 VELOCITY AND SLOPE FACTOBS 

There are several ultrasonic attenuation 
factors that correlate with naterial 
properties. These factors derive froa 
the AFC paraaeters (C) and (B) , as 
explained next. 

Since attenuation is a strong function 
of frequency, it is necessary to specify 
particular frequencies that are relevant 
for aaterial property correlations. 
There are two factors based on the first 
derivative or slope of the ifC at 
specific frequencies: One factor is the 
slope (BETA1) which is d (ALPHA) /d (F) 
evaluated at the frequency corresponding 
tc an attenuation coefficient of unity 
or, (ALPHA) - 1. The other factor is the 
slope (BETA2) which is d (ALPHA) /d (F) 
evaluated at a frequency corresponding 
to a particular wavelength. This is the 
wavelength aarking the transition froa 
Rayleigh to stochastic scattering. It 
is defined in teras of the Bean grain 
size (GS) of the aaterial. Expressions 
for these factors are<*>, 

LOGN (BETA1) * LOGN(M) ♦ [ LOON(C)] /(M) 

LOGN (BET A2) = LOGN[ (H) (C) ] * 

[ (H)-1](L0GN[ (VL)/(GS) ]) 

According to refs. 2,5 the quantity 
(BETA 1 ) increases linearly with plane 
strain fracture toughness (KIC) for 
aaterial saaples of a aaraging steel if 
all saaples have the sane yield strength 
(YS). The product of velocity with 
(BET42)/ (H) , that is, (VL> (BETA2) / (H) , 
varies exponentially with the quantity 
(KIC) */ (IS) * for the aaraging steel. 
These relations are shown in figs. 12(a) 
and 12(b). In fig. 12 data obtained by 
digital wavefora processing are conpared 
with previous data obtained by analog 
aethods described in ref. 2. 

5.4 DISCUSSION 

Ultimately, the usefulness of any set of 
velocity or attenuation aeasureaents or 
derived quantities aust be judged on the 
basis of capability for characterizing 
aaterial properties with good accuracy. 


The digital wavefora processing aethods 
described herein appear viable froa this 
standpoint. For exaaple, the data in 
figs. 12(a) and 12(b) confira the strong 
correlations that were expected between 
the ultrasonic factors and fracture 
toughness and yield strength. 

In figs. 12(a) and 12(b) the slight 
difference between the regression lines 
for digital and analog data is due 
priaarily to differences in the 
reflection coefficients used. In the 
aanual analog netbod the reflection 
coefficient was taken as (BCC) , because 
it was inferred to be appropriate based 
on the available data. The autoaatic 
digital aethod, on the other hand, 
produces considerably aore data and 
better resolution in the frequency 
doaain. It becoaes clear that (BCC) is 
only an upper bound on the reflection 
coefficient. 

It is noteworthy that the ARC intercepts 
the ordinate asyaptotically to (BCV) 
between (BCC) and (RCU) , as in figs. 8 
and 10. (BCV) which is based on the ABC 
is thus the "valid" reflection 
coefficient relative to the high 
frequency data in the valid zone. 
Horeover, when good coupling is attained 
by applying adequate force, (RCV) tends 
closer to (RCB) . This is expected if 
the buffer is tightly coupled to the 
aaterial sample surface. 

Havefora truncation results revealed an 
interesting aspect of the trailing 
oscillations in the natural wavefora. 
The truncation criterion was to retain 
only the five syaaetrical leading 
oscillations, and to discard the snail 
trailing oscillations. Once a wavefora 
aeets the acceptance criteria, given 
previously, it is a simple natter to 
perform the truncation precisely and 
reproducibly using digital wavefora 

processing. The discarded trailing 

oscillations of the natural wavefora are 
assuaed to be due chiefly to nultiple 

residual reflections in the couplant. 
This is confiraed by comparing the ABC 
raw data for natural and trucated 
waveforms as follows: 

The raw data in figs. 10(a) and 10(b) 
reveal the effects of truncation in the 
diffraction regine, to the left of the 
valid zone. it is apparent in fig. 

10 (b) that truncation eliainates the 
effects of couplant nultiple reflection 
seen in fig. 10 (a). 

These multiple reflection effects are 
apparently due to couplant thickness 
variations since the surfaces are 
neither prefectly sated nor ideally 
planar. Hence, in the low frequency 
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revise, the loncer wavelengths are 
likely to interact with a greater 
portion of the coupled area. The result 
is the anoealous behavior of the raw 
data in the diffraction regiae. In the 
diffraction regiae the reflection 
coefficient (BC) is a variable that 
depends on the coupling nediun, 
wavelength, and ultrasonic bean cross 
section* 7 >. The anonolous reflection 
effects can be avoided as denonstrated 
by the ABC which is based only on data 
in the high frequency regiae (valid 
zone). The reflection coefficient (BCV) 
is thus established independently of the 
anoaalous effects in the diffraction 
regiae. 

Clearly, autonated digital waveforn 
processing nethods will provide the 
experiaentalist additional inforaation 
that is either unavailable or 
iapractical to obtain with aanual analog 
nethods. This additional inforaation is 
gained through the coaputer's ability to 
process aassive aaounts of data. This 
affords access to data that is readily 
available only by aeans of conplex 
analytical treataents such as the fast 
Fourier transfora and array processing 
routines. 

6. summary and conclusions 

1 coaputerized waveforn processing 
systea and its operational details are 
described. The systea was dedicated to 
ultrasonic velocity and attenuation 
aeasureaents relative to evaluation of 
aaterial properties. The systea serves 
sinultaneously as a tool for ultrasonic 
research and for establishing aaterial 
property correlations. 

The systea and nethods described herein 
are essentially prototypical. They 
denonstrate soae practical approaches to 
ultrasonic nondestructive aaterials 
evaluation. More work is needed to 
iaprove the operating processes and 
instruaentation. Aaong the advancenents 
needed are: 

(1) Faster algo rl thas in nachine 
language for soae aspects of data 
processing not currently covered 
by existing software 

(2) Fully progranaable iastruaents 
and nodules for ultrasonic signal 
generation, acquisition, and 
aaplif icatioa 

(3) Bore powerful and nore sensitive 
transducers and instruaents for use 
over greater frequency ranges 

(to several hundred negahertz) 


It is apparent that high speed digital 
waveforn processing over a wider 
frequency range will enable lore 
detailed ultrasonic aeasureaents for 
aaterial property evaluations. Siailar 
aeasureaent have been previously 
accoaplished by aanual analog aethods 
that are tedious and tine consuaing 
while giving less detailed data. 
Without coaputer autoaation the 
necessary aeasureaents for aaterial 
property evaluations would certainly be 
iapractical either in production or 
laboratory situations requiring speed, 
ecoaoay, and accuracy. 

The experience gained thus far indicates 
that further developaent of coaputer 
wavefora processing in high frequency 
ultrasonics will contrioute to: 

(1) Establishing a powerful research 
tool for ultrasonic investigations. 

(2) Practical applications for aaterial 
property verification and control. 

(3) Investigation of factors governing 
aaterial strength. 
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SPECIMEN' 33-2 
IE-3 U 



(a) ECHO OVERLAP OF B1 AND B2. 



FREQUENCY, MHZ 

lb) SPECTRA FOR B1 AND B2. 


PRC PARAMETERS 
ZONE' 20 - 37 NH2 
<M>- 3 03 
<C>- 1 3940 2E -03 
( FIT >• .998929 



FREQUENCY, tWZ 

ic) LINEAR PLOT OF (ALPHA) VS. |F). 


Figure 11. - Computer documentation of data for maraging steel specimen 33-2. Result 
of echo overlap routine used for computing velocity is shown in (a). Valid zone for 
determining ARC parameters is indicated by vertical markers in (c). Final documen- 
tation of input data, conditions, and processing results is illustrated in Id). 


SPECIMEN 

MATERIAL 


33-2 

MS-230 


ATTENUATION US. FREQUENCY CURUE 


THICKNESS < ST )■ 2492 CM 

DENSITY < ON >» 8 03 G^CC 
GRAIN SIZE <GS>« 13 UM 

UELOCITY < UL >“ 367826 CIUUS 

CENTER FREQUENCY (CF>- 34 MHZ 

TRUNCATE!) 

ARC PARAMETERS' 

UAL ID ZONE 20 TO 37 MHZ 
REFLECTION <RCU» 396138 
<M>- 3 05 
<C» 1 39402E-83 
<FIT>- 998929 

*FC PARAMETERS' 

< BETA1 >“ 0780292 

( BETA2 >• 10 9933 
<ULXBETA2>- 6 24226 
<ULXBETA2V<M>- 2 04664 

TRANSDUCER 30 MHZ 
TYPE' LONGITUOINAL 
COUPLANT' GLYCERINE 

REPETITIONS^POSITIONS- 3 ✓ 1 



Id) SPECIMEN DATA AND LOGARITHMIC PLOT OF (ALPHA) VS. (F). 


Figure 11. - Concluded. 
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ial PLOT OF (BETAll VALUES VS. FRACTURE TOUGHNESS. 



(bl PLOT Of (VLMBCTA2I iM» VALUES VS. CHARACURISTIC LENGTH. 


figure 12. Comparison of ultrasonic attenuation factors with 
theoretical curves (solid lines > tor tour maragmg steel 1250- 
grade* specimens. 



